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ABSTRACT 

We present interferometric imaging at 33 GHz of the Corona Borealis supercluster, us- 
ing the extended configuration of the Very Small Array. A total area of 24 deg 2 has been 
imaged, with an angular resolution of 11 arcmin and a sensitivity of 12 mJy/beam. 
The aim of these observations is to search for Sunyaev-Zel'dovich (SZ) detections 
from known clusters of galaxies in this supercluster and for a possible extended SZ 
decrement due to diffuse warm/hot gas in the intercluster medium. Hydrodynamical 
simulations suggest that a significant part of the missing baryons in the local Universe 
may be located in superclusters. 

The maps constructed from these observations have a significant contribution 
from primordial fluctuations. We measure negative flux values in the positions of the 
ten richest clusters in the region. Collectively, this implies a 3.0-sigma detection of 
the SZ effect. For two of these clusters, A2061 and A2065, we find decrements of 
approximately 2-sigma each. 

Our main result is the detection of two strong and resolved negative features at 
-70 ± 12 mJy/beam (-157 ± 27 /zK) and -103 ± 10 mJy/beam (-230 ± 23 /kK), 
respectively, located in a region with no known clusters, near the centre of the su- 
percluster. We discuss their possible origins in terms of primordial CMB anisotropies 
and/or SZ signals related to either unknown clusters or to a diffuse extended warm/hot 
gas distribution. Our analyses have revealed that a primordial CMB fluctuation is a 
plausible explanation for the weaker feature (probability of 37.82%). For the stronger 
one, neither primordial CMB (probability of 0.33%) nor SZ can account alone for its 
size and total intensity. The most reasonable explanation, then, is a combination of 
both primordial CMB and SZ signal. Finally, we explore what characteristics would 
be required for a filamentary structure consisting of warm/hot diffuse gas in order to 
produce a significant contribution to such a spot taking into account the constraints 
set by X-ray data. 

Key words: cosmology: observations - cosmic microwave background - galax- 
iesxlusters - techniques:interferometric 
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1 INTRODUCTION 

The mean baryon density of the Universe is one of the most 
relevant cosmological parameters, as it influences baryonic 
structures on all scales, from the abundances of primordial 
nuclei to the large-scale distribution of galaxies and inter- 
galactic gas. Indeed, its value is also a prediction of the stan- 
dard Big Bang model that may be tested observationally. 
The calculation of the primaeval abundances of light ele- 
ments, such as deuterium, combined with standard Big Bang 
nucleosynthesis, allows a very pr ecise determination of this 
param eter for standard models llBurles. Nollett. fc Turner! 
l200ll) : n B = (0.020±0.002) h~ 2 = 0.038±0.004, where here 
h = ftioo = .Ho/100 km s~ 1 Mpc~ 1 = .73 is adopted for the 
last term. Likewise iRauch et alJ (Il997f) derive a lower limit 
fl B = 0.021 h~ 2 = 0.042 for a ACDM model from obser- 
vations of the Lya forest absorption in a selected sample of 
seven high resolution quasar spectra at z = 2. More recently, 
the best-fit cosmological model to the first year WMAP data 
release also gi ves a value n B = ( 0.0224 ± 0.0009) h~ 2 = 
0.044 ± 0.002 JSpergel et alJ 120031) . while the recent VSA 
results give Q B = (0.0234±g;ggi|) h~ 2 = 0.0464±°S 
jRebolo et alJl2004h . 

In principle, the consistency between these three com- 
pletely independent methods of estimating the baryon 
density is straightforward, but at z = in the 
present day Universe, the sum over all the well observed 
components give a considerably smaller value. Indeed, 
iFukueita. Hoean. fc Peebles! <ll998T) made an estimate of the 
global budget of baryons in all known states, namely the 
different kinds of stars, the atomic and molecular gas, and 
the plasma in galaxy clusters and in groups. They infer a 
value Q B = (0.010 ± 0.003) h~ 2 = 0.020 ± 0.007. 

Thus, if no serious errors have been made in the theo- 
retical estimates, which seems unlikely given the agreement 
between the different methods, it seems that most of the 
baryons are yet to be detected in the present day Universe. 
This is the well-known "baryon problem" . 

One of the most important hypotheses related to this 
hidden matter is based upon hydrodyn amical simulations 
<Cen fc OstrikeJ Il999t iDave et al.1 l200ll) that predict the 
formation at low redshift (z < 1) of a very diffuse gas 
phase with temperatures 10 5 < T < 10 7 K, which is 
neither low enough to have permitted condensations into 
stars or to form cool galactic gas, nor as high as that of 
the hot gas present in galaxy clusters. Such low density 
gas could account for a substantial fraction of the missing 
baryons, and according to these simulations it should be 
distributed in large scale sheet-like structures and filaments 
connecting clusters of galaxies, due to the infall of bary- 
onic matter into previously formed dark matter filaments. 
Moreover, this amount of gas would not violate the con- 
straints on the spectral distortions of the cosmic microwave 
background (CMB). It constitutes what is known as the 
"warm/hot intergalactic medium" (WHIM) and would be 
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observed in the soft X-ray band. The detection of its ra- 
diation could be obscured by the presence of many galac- 
tic foregrounds and by extragalactic contributions from 
groups of galaxies, clusters or AGNs. Nevertheless, sev- 
eral attempts in recent years have been made to detect it, 
and some detections have been claimed, either by studying 
the correlation between the observed so ft X-ray structures 
and selected galaxy ov erdense regions dScharf et alJ 12000: 
IZappacosta et alJl2002[) or bv detecting a soft X-ray ex cess 
in clusters of galax ies iFinoguenov. Briel. fc Henrvteooljl) . or 
in their proximit y jBriel^^H^rirv^^9^l TrttI^v^^H^ririksenl 
l200lHSoltan. Frevberg fc HasingeiH2002f) . 

As indicated by ICen fc Ostrikej {l999'), both X-ray 
emission and the thermal Sunyaev-Zel'dovich (SZ) effect 
can be cross-correlated with galaxy or galaxy cluster 
catalogues to search for the missing baryons. The SZ effect 
typically originates from the richest clusters of galaxies that 
contain extended atmospheres of hot gas (k B T e ~ 10 keV). 
But there may be other objects that could also produce 
a detectable SZ signal, such as superclusters of galaxies, 
where even though low baryon overdensities are expected, 
path lengths may be long so that a significant SZ effect 
could build up since the effect is propo rtional to the lin e 
of sight integral of the electron density feirkinshawlll999T) . 
Also, it is reasonable to look at superclusters as optimal 
regions for the detection of WHIM, since simulations 
show that this gas should be distributed in filamentary 
structures extending over several t ens of Mpc and con- 
nectin g clusters of galaxies. In fact, iMolnar fc Birkinshawl 
(1998) looked in the COBE-DMR data in the region of 
the Shapley supercluster, but found no clear evidence of 
an extended SZ effect and therefore set an upper limit of 
< —100 /iK on the gas temperature. On the other hand, 
several studies have been carried out not by searching 
in selected overdense regions, but by trying to obtain a 
statistical SZ detection from intra-supercluster (ISC) gas 
from analyses of the whole sky, first using the COBE-DMK 
data feandav et alJ Il996[) and, more recently, the first 
relea se of WMAP data ^osa^ba^Ga^jsSMS^-^LS^ i^idei 
1 20031: Illernandez-Monteagudo fc Rubino-Martml 120041 : 
iMvers et alJ 120041 : iHernandez-Monteagudo et alJ l2004h . 
This is mainly carried out by correlating CMB maps with 
templates built from galaxy or galaxy cluster catalogues. 

In this study we concentrate on the Corona Borealis 
supercluster (CrB-SC). This supercluster does not contain 
high flux density radio sources which may inhibit SZ mea- 
surements. T he supercluster members have relatively high 
X-ray fluxes dEbeling et aljll998l l2000F) which suggests that 
detectable SZ effects could build up. Also, the angular size 
of the supercluster core is suitable for mosaiced observations 
with the VSA extended array, which has a primary beam 
FWHM of 2.1°. Taking these factors into account, we se- 
lected this supercluster as the most interesting object of its 
type within the region of sky observable by the VSA (which 
have a declination coverage of ~ — 7° to +63°). 

In Section 2 we present a general description of the CrB- 
SC and summarize previous observational work. Section 3 
includes an overview of the VSA interferometer and explains 
the data reduction and the map making procedure. In Sec- 
tion 4 we discuss the SZ effects from the individual clusters 
in the region and the possible causes of the two observed 
negative features. Conclusions are presented in Section 5. 
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Figure 1. Spatial distribution of the Abell and Zwicky clusters 
in the CrB-SC core region. The Abell clusters labelled with blue 
letters belong to the CrB-SC. Also indicated are the positions 
of the S CL168 (CrB) and SC L157 supercluster centres, as re- 
ported bv lEinasto et all Red dots indicate the positions of 
the galaxies at redshifts within ±0.02 of the supercluster redshift 
and inside the optical radius of each Abell cluster. (The redshifts 
and the optical radius are taken from the NASA Extragalactic 
Database.) 



study of a sample of 1555 galaxies in the vicinity of Abell 
clusters. They report 97 new redshift measurements over the 
previous 85. They conclude that the masses of all clusters in 
the core of the CrB-SC lie in the range 1.5-8.9 x 10 14 M , 
while the mass of the SC is approximately « 8.2 x 10 Mq, 
which is probably enough to bind the system. They proved 
also that the dynamical timescales are comparable with 
the Hubble time, making it unlikely that the system could 
be virialis ed, as one might expect. This work was ex- 
tended by ISmall. Ma. Sargent, fc Hamilton! <199Sft by in- 
creasing the number of galaxies with known redshifts to 
528. They quote a value for the mass of the supercluster 
of 3 - 8 x 10 16 h' 1 Mq, a v alue slightly higher than the 
previo us one. This is because IPostman. Geller. fc Huchral 
( 1988) made the assumption that the mass-to-light ratio 
in the CrB-SC is the same as in the richest Abell clusters 
and they used a supercluster volume three times smaller. 
They remark that almost one third of the galaxies in the 
region are not linked to any Abell cluster, and also empha- 
size the great contribution to the projected surface density 
of galaxies of the background cluster A2069 and its sur- 
rounding galaxies, located at a redshift z « 0.11, suggest- 
ing the existenc e of the so-cal l ed "A 2069 supercluster" . On 
the other hand, iMarini et alJ i2004l) have analysed on the 
BeppoSax X-ray data in a search for evidence of a merging 
signature between the pairs of clusters candidates in the re- 
gion of CrB-SC A2061-A2067 and A2122-A2124. They find 
no clear evidence of interaction but detect a candidate shock 
inside A2061. 



2 THE CORONA BOREALIS SUPERCLUSTER 

The CrB-SC is one of the mos t prominent examples of su per- 
clustering in the northern skv. lShane fc Wirtanenl l|l954f) . af- 
ter counting galaxies on the Lick Observatory photographic 
plates, were the first to remark upon the extraordinary cloud 
of gala xies that constitute the supercluster, and later lAbelll 
( 1958) noted the presence of this concentration of clusters of 
galaxies including it in his catalogue of second order clusters. 
Depending on author, the num ber of clusters belonging to 
this SC ranges from six to eight ([ Postman" Geller. fc Huchral 
Il988l: ISmall Sargent, fc Hamilton: ! 11997]), but w e will focus 
on the classification given in the lEinasto et alJ <|200 ll) cat- 
alogue, according to which the CrB-SC includes eight clus- 
ters, around the position R.A. = 15 h 25 m 16.2 s , DEC = 
+29°31'30", at a redshift z « 0.07. These clusters are listed 
in Table along with their characteristics. Six of these 
(A2061, A2065, A2067, A2079, A2089 and A2092) are lo- 
cated in the core of the SC, in a w 3° x 3° region, while 
there are two others (A2019 a nd A2124) at an ang ular dis- 
tance of w 2.5° from the core. lEinasto et al.l B2001I) include 
only A2061 and A2065 as X-ray emitting clusters. There are 
four other Abell clusters (A2056, A2005, A2022 and A2122) 
in this region at redshifts around z fa 0.07, two (A2069 and 
A2083) at z « 0.11 and another two (A2059 and A2073) 
that are even more distant. There are also several Zwicky 
clusters but without redshift measurements. The spatial dis- 
tribution of these clusters can be seen in Figure 

The fi rst dynamical study of the Cr B-SC was car- 
ried out bv lPostman. Geller. fc Huchral <ll988l) . through the 



3 OBSERVATIONS AND DATA REDUCTION 

3.1 Observations with the VSA interferometer 

The Very Small Array (VSA) is a 14-element heterodyne 
interferometric array, tunable between 26 and 36 GHz with 
a 1.5 GHz bandwidth and a system temperature of approx- 
imately 35 K. It is located at an altitude of 2400 m at the 
Teide Observatory in Tenerife. For this study, the observ- 
ing frequency was set at 33 GHz, and we used the extended 
configuration, which uses conical corrugated horn antennae 
with 322 mm apertures, and has a primary beam FWHM 
of 2.1°, and a synthesised beam FWHM w 11 arcmin at 
the observing freq uency. For a detailed description of the 
instrument see e.g. lWatson et all (|2003). 

Situated next to the main array is a two-element in- 
terferometer, which consists of two 3.7-m diameter dishes 
with a north-south baseline of 9 m, giving a resolution of 

4 arcmin in a 9 arcmin field. This interferometer is used for 
source subtraction (SS), monitoring radio sources simulta- 
neously with the m ain array observations, as described in 
iTavlor et alJ (|2(]03). This strategy neatly copes with the ad- 
ditional problem of source variability. 

The observations were carried out during the period 
spring 2003 - summer 2004, although most of the data were 
taken in the first months, while the last months were used 
only to improve the noise levels in those fields with a poorer 
signal to noise, and to obtain a better measurement of the 
fluxes of some of the radio sources. Initially, we observed five 
individual patches covering the interesting central regions of 
the supercluster, plus two additional pointings separate from 
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Table 1. Clusters: basic data. The coordinates are the optical centre of each cluster. The last two columns show the predicted flux 
which would be observed by the VSA for those clusters with X-ray measurements (labelled "pred", see Section 4.1 for details), and the 
measured value in our final maps (labelled "VSA"). The error bars contain the thermal noise (<r n ), the average primordial CMB (which 
has been computed to be ocmb ~ 19 mjy/beam) and the residual sources (cr sour rj 6 mjy/beam) added in quadrature. 



Cluster 


RAM 


DECM 


zO) 


£40.1-2.4 keV) 




AS pred 


AS^ SA 




(J2000) 


(J2000) 




(10 44 h~ 2 ergs- 1 ) 


(keV) 


(mjy/beam) 


(mjy/beam) 


A2019 


15 02 57.2 


+27 11 17 


0.0807 








-2.6 ± 24.5 


A2061 


15 21 15.3 


+30 39 17 


0.0784 


3.950) 


5.6( b ) 


-54.7 


-53.8 + 27.4 


A2065 


15 22 42.6 


+27 43 21 


0.0726 


4.94(6) 


8.4W 


-61.7 


-42.8 + 23.4 


A2067 


15 23 14.8 


+30 54 23 


0.0739 


0.86( c ' 


3.1< c ' 


-27.7 


-26.5 ± 27.4 


SCL-158 


15 25 16.2 


+29 31 30 












A2079 


15 28 04.7 


+28 52 40 


0.0690 








-6.7 + 22.6 


A2089 


15 32 41.3 


+28 00 56 


0.0731 








-28.9 + 27.9 


A2092 


15 33 19.4 


+31 08 58 


0.0669 








-21.0 + 28.5 


A2124 


15 44 59.3 


+36 03 40 


0.0656 


1.35W 


3.7W 


-44.3 


-19.4 + 29.2 


A2069 


15 23 57.9 


+29 53 26 


0.1160 


8.92<» 


7.9< 6 > 


-57.1 


-15.2 ± 27.4 


A2073 


15 25 41.5 


+28 24 32 


0.1717 


3.74( c ) 


5.6^ 


-29.9 


-33.4 + 22.6 



Notes: (a) NASA Extragalactic Database, Cbl lEbeling et al] <1998L Ccl lEbeling et al.l i2000l) . 

Note the use of /150 = Hg/50 km s — 1 Afpc - 1 . 



the main mosaic and centred on clusters A2019 and A2124. 
We later added four more pointings to improve the coverage 
of the supercluster core. These observations are described 
in Table [5] Each daily observation was approximately four 
hours in duration. We dedicated between 5 (for CrB-D and 
CrB-G) and 37 (for CrB-H) days to each pointing. We quote 
both total observation and integration times, the latter indi- 
cating the data retained after flagging (for instance, periods 
of bad weather). 



3.2 Calibration and data reduction 

The primary calibrator for the VSA is Jupiter. The flux 
scale is transferred to other calibration sources: Tau A (the 
Crab Nebula) and Cas A, which are observed daily. A de- 
tai led description of the VSA calibrati on process is pr e sented 
in I Watson et all (120031) . w hereas in iGrainge etafl (120031) 
and TPickinso'ii'et al ] J2004I) are remarked the specifications 
adopted for the VSA extended config uration. 

In the first VS A studies (e.g. IWatson etaD l)2003ft : 
IGrainge et al.l teOOStO the calibration scale was based on 
the ef fective temperature of Jupiter given by iMason et al] 
( 1999), which extrapolated to 33 GHz is Tj up = 153 ± 5 K 
(3% accurac y in temperature). In the most recent VSA pub- 
lished data ^Dickinson et al.ll2004h we rescale our results us- 
ing the Jupiter temperature as derived from the WMAP 
first- year data: T Jup = 146.6 ± 2 K JPage et al-lEooll) . 
which agrees with the lMason et, al.l l|l999T) value at la level, 
but reduces the calibration error from 3 to 1.5% in tempera- 
ture terms (6% to 3% in the power spectrum). In this paper 
we have also rescaled the Jupiter temperature to the WMAP 
value, which is the most accurate value published to date. 

The correlated signal from each of 91 VSA base- 
lines is proces s ed as described in IWatson et al.l (I2003T) and 
IGrainge et al.l (120031) . and the data checks so described are 
also applied. 



3.3 Source subtraction 

In order to remove the effect of radio sources from the 
data, we followed a slightly different approach to that con- 
sidere d for the primordial CMB fields (see IClearv et al] 
J2004I) '). For those fields, the Ryle Telescope (RT) was 
used to iden tify all the sources above a given threshold at 
15 GHz (see iTavlor et all feOOSt and IGrainge et all (120031) 
for the compact and extended arrays, respectively). Us- 
ing this catalogue, the source subtraction baseline moni- 
tored all these sources in real time, and the derived fluxes 
were subtracted from the VSA main array data. How- 
ever, for the present study a complete scanning of the re- 
gion with the RT was not available owing to observing 
time constraints, so in order to build a source catalogue 
in the region, we proce eded as follows. All the sources 
in the N VSS-1.4 GHz JCondon et al] Il998l) and GB6- 
4.85 GHz dGregorv. Scott. Douglas, fc Condonlll99rj) cata- 
logues closer than 2° from each pointing position were iden- 
tified. Using the fluxes from these catalogues, the inferred 
spectral index between the respective frequencies was used 
to extrapolate the flux to 33 GHz. Finally, all sources with 
predicted fluxes above 20 mjy at 33 GHz were monitored 
using the SS baseline. The flux threshold was chosen accord- 
ing to the requirements for the primordial CMB fields, but, 
as we shall see below, for the purposes of this paper this 
limit can be relaxed. 

Source positions were assigned according to the coor- 
dinates in the GB6 catalogue. This catalogue has a lower 
angular resolution (FWHM ~ 3.5 arcmin) than the NVSS 
(FWHM ~ 45 arcsec), so multiple sources in NVSS could 
be associated with a single entry in the GB6 catalogue. This 
was taken into account when computing the spectral index 
for each source: the adopted flux at the NVSS frequency was 
derived using the GB6 resolution, combining all sources in 
NVSS which correspond to a given entry in the GB6 cat- 
alogue (i.e. for each GB6 source, we found all the NVSS 
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Table 2. Summary of the observations with the VSA. We list the central coordinate of the eleven pointings, the 
total observation time, the integration time and the thermal noise achieved for each map (computed from the map 
far away from the primary beam). 



Pointing RA (J2000) DEC (J2000) T obs T int Thermal noise 

(hr) (hr) (mjy/beam) 



CrB-A 


15 23 12.00 


+28 06 00.0 


54 


50 


12.4 


CrB-B 


15 27 48.00 


+29 24 00.0 


70 


70 


10.8 


CrB-C 


15 22 48.00 


+30 21 00.0 


33 


33 


18.9 


CrB-D 


15 32 00.00 


+30 45 00.0 


19 


19 


20.5 


CrB-E 


15 32 00.00 


+28 18 00.0 


22 


22 


19.7 


CrB-F 


15 02 57.20 


+27 11 17.3 


47 


43 


14.4 


CrB-G 


15 45 00.00 


+36 03 57.6 


19 


19 


21.4 


CrB-H 


15 23 00.00 


+29 13 30.0 


167 


130 


10.2 


CrB-I 


15 27 24.00 


+30 33 00.0 


56 


11 


18.6 


CrB-J 


15 32 00.00 


+29 31 30.0 


55 


39 


18.9 


CrB-K 


15 28 00.00 


+28 12 00.0 


41 


33 


20.3 



Table 3. List of all the radio sources with measured fluxes at 33 GHz above 50 mjy. We quote the extrapolated 
flux at 33 GHz, and the estimated flux at 33 GHz using the SS observations, which have been used to subtract the 
sources from the data. 





Name 


RA (J2000) 


DEC (J2000) 


Extrapolated flux 


Measured flux by 










NVSS-GB6 at 33 GHz (mjy) 


the SS at 33 GHz (mjy) 


CrB-F 


1459+2708 


14 59 39.60 


27 8 16.0 


83 


132 ± 12 




1504+2854 


15 4 27.30 


28 54 25.0 


76 


193 ± 14 




1509+2642 


15 9 39.10 


26 42 45.0 


39 


58 ± 10 


Mosaic 


1514+2931 


15 14 20.90 


29 31 9.0 


104 


95 + 6 




1514+2855 


15 14 40.30 


28 55 39.0 


72 


63 + 8 




1514+2943 


15 14 3.70 


29 43 21.0 


33 


50 ± 13 




1521+3115 


15 21 1.80 


31 15 50.0 


196 


69 ±43 




1522+3144 


15 22 9.50 


31 44 18.0 


243 


296 ± 7 




1522+2808 


15 22 48.90 


28 8 51.0 


133 


97 + 5 




1527+3115 


15 27 18.20 


31 15 14.0 


294 


157 + 6 




1528+3157 


15 28 52.60 


31 57 34.0 


75 


98 + 7 




1529+3225 


15 29 38.70 


32 25 23.0 


60 


84 ±38 




1531+2819 


15 31 21.40 


28 19 26.0 


30 


58 ±9 




1532+2919 


15 32 20.20 


29 19 40.0 


71 


51 ± 12 




1535+3126 


15 35 58.90 


31 26 25.0 


47 


53 ±5 




1537+2648 


15 37 6.40 


26 48 24.0 


23 


56 ± 18 




1539+3103 


15 39 15.90 


31 3 59.0 


107 


86 ± 15 




1539+2744 


15 39 38.80 


27 44 33.0 


244 


218 ±10 


CrB-G 


1538+3557 


15 38 57.40 


35 57 9.0 


37 


54 ± 12 




1540+3538 


15 40 31.70 


35 38 26.0 


24 


52 ± 16 




1544+3713 


15 44 44.70 


37 13 22.0 


25 


50 ± 11 




1546+3631 


15 46 7.60 


36 31 7.0 


26 


51 ± 18 




1546+3644 


15 46 38.30 


36 44 30.0 


44 


59 ± 15 




1547+3518 


15 47 51.70 


35 18 59.0 


87 


77 ± 12 




1552+3716 


15 52 5.30 


37 16 5.0 


47 


50 ± 18 



sources closer than 5 arcmin, applied the primary beam cor- 
rection, and summed the corrected fluxes). 

Taking into account that the sensitivity of the GB6 
catalogue is 18 mjy at 4.85 GHz, together with a maxi- 
mum value o f the rising spectral i ndex between 1.4 GHz and 
33 GHz (e.g. iMason et all Jl999M of a = 0.5 (S oc u a ), we 
are confident that we have measured all radio sources with 
fluxes above ~ 50 mjy at 33 GHz. Nevertheless, we appreci- 
ate that sources with inverted spectra and higher fluxes may 



exist, and identify four compact features at the 3a n level (be- 
ing a n the standard deviation of the thermal noise) in our 
maps of the CrB-SC region. We suspected that these may 
be radio sources not picked up by our strategy, but source 
subtractor observations suggest otherwise. No sources with 
detectable fluxes (i.e. S > 20 mjy at the 2a level) were 
found at those positions. 

A total of 74 sources closer than 2° to any pointing were 
found to have an extrapolated flux greater than 20 mjy at 
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33 GHz, and of these, only 25 had a measured flux with 
the source subtractor greater than 50 mjy. In Table [3] wc 
present the final list of identified sources with measured 
fluxes greater than 50 mjy at 33 GHz. These values were 
used to carry out the source subtraction. 

It should be noted that this flux limit is larger than 
that used for the analyses of the primordial CMB fields 
with the extended VSA array. For that study, we used 
a flux limit of 20 mjy, which guarantees that the resid- 
ual contribution from unsubtracted s ources is less than the 
flux sensitivity l|Grainge et aljl2003l : iDickinson et al.ll2004t 
IClearv et alll2004) . and this could be achieved thanks to 
the RT survey at 15 GHz. However, as we shall see below 
(Section 4.2), the main result of this paper is the detection 
of two anomalous cold regions with the VSA in the CrB-SC 
region, that seem difficult to explain as primordial CMB fea- 
tures. We took our list of extrapolated fluxes and used the 
values to determine the effects on our measurements of the 
cold spots. We found that the individual effect of any one 
source was sufficient only to affect our measurements by less 
than 2%, and that the collective effect of all sources would 
cause effects of 1.3% and 3.2% in the two spots respectively. 
This latter error is only twice the calibration error. Thus 
for the purposes of this paper (in contrast with the VSA 
power spectrum measurem ents, where more stringent source 
constraints were required llScott et al.ll2003t iGrainge et alJ 
120031: IDickinson et al J 12004) 1. our source subtraction limit 
of 50 mjy is sufficient. 

In any case, an estimate of the confusion noise due to 
unresolved sources below our source subtraction threshold 
can be made using the best fit power-law model obtained 
from the source count above ^ 50 mjy (see fluxes in Ta- 
ble EJ, 



n(S) = (17.6±2.3) 



S 



70 mjy 



mjy 1 sr 1 . (1) 



Note that this fit has a slope comp atible with that of the 
model presented in lClearv et alJ <l2004l) . which was obtained 
using the source count from the VSA primordial CMB obser- 
vations. However, the amplitude in our case is a factor w 1.7 
larger. This is as expected since there may be a higher den- 
sity of radio sources in the supercluster due to the presence 
of source s associated with th e member clusters. Indeed, the 
model of IClearv et al] <l2004f) predicts just 15 radio sources 
above 50 mjy in our region of observation, whereas we find 
25. Using the model described by equation we estimate 
that the residual sources may introduce a confusion noise 
o~sour = 7.3 mjy/beam. Note however that this result may 
overestimate the real confusion noise in the regions far away 
from the CrB-SC clusters, where the source count sh ould be 
described by the model given in IClearv et al] ]2004l) . 

To conclude this section, we also report a problem which 
was detected during a late stage in the data processing, re- 
lated to SS measurements. Since the source subtractor con- 
tinually slews between pointed observations of point sources 
inside the observed VSA fields, the phase stability of the 
system has to be che cked several times du ring an observing 
run, as described in IWatson et al.l <l2003l) . To this end, we 
observe several interleaved strong calibration sources close 
to the considered region on the sky in order to track the 
phase stability of the system properly. We found that one 



of the calibration sources used in the CrB-SC central mo- 
saic was not as bright as predicted by the extrapolation and 
hence was not bright enough to produce a sufficiently ac- 
curate determination of the phase. Using Monte Carlo sim- 
ulations with the known calibrator flux and the sensitivity 
of the measurement, it is easy to show that a weak calibra- 
tion source produces an average underestimate of the source 
flux. From these simulations, we could try to correct all the 
flux measurements, but this would require us to make an 
assumption about the variability of each source. Instead, we 
decided to use a new estimator of the flux of the source that 
does not require phase calibration. 

Let us introduce the estimator E[S 2 ] — Vr + V 2 — 2a 2 
where S is the flux of the source we are measuring, Vr and Vi 
the real and the imaginary parts of the visibility, and a the 
statistical noise of each measurement. By construction, this 
estimator is unbiased ({E[S 2 ]} = {S 2 }) and is not sensitive 
to phase error of the instrument. If the noise is Gaussian, 
then E is distributed as a (displaced) \ 2 with two degrees 
of freedom. However, we must note that for our purposes, 
this estimator can only be used for the case of non-variable 
sources, because it provides (S 2 ) during the whole period 
of observation, and not (S), which is the quantity we are 
interested in. Thus, we finally decided to use the estimator 
E' — vE to derive the fluxes. This estimator is obviously 
defined only for positive values of E and in principle is bi- 
ased to high values of the flux if the signal-to-noise ratio 
of the measurement is small. This bias was computed and 
we find that, for the flux sensitivity of the instrument, it is 
significant only for values of the flux smaller than 50 mjy. 
However, we have reduced this bias by averaging the square 
estimator, E, between adjacent measurements of the same 
source to get a positive value prior to application of the 
square root. Summarizing, even if we do not have an accu- 
rate determination of the phase of the instrument at the time 
of the observation for some of the sources, their fluxes can 
be recovered from our measurements using this new estima- 
tor, although the error bar on the final measurement will be 
larger than the one obtained with the standard procedure. 

In addition, and as a consistency check, we re-observed 
the strongest sources in the mosaic, with estimated fluxes 
above w 60 mjy, but now using an appropriate calibrator. 
A variability study was then ap plied to all t h ese so urces 
following the same criteria as in IClearv et al] J2004), and 
variable sources were identified. Only five sources were found 
to be variable. For non-variable sources, the new flux was 
directly compared to our previous determination to check 
that our flux values are robust. Final values for the fluxes 
are presented in Tabled 



3.4 VSA maps 

Daily observations are calibrated and reduced as described 
in Section 3.2 (and references therein) and are held as visi- 
bility files, which contain the real and imaginary parts mea- 
sured at each uu-position together with the associated rms 
noise level. T hese files are loaded individually into AIPS 
iGreisenlll994l) . where the map making process is carried 
out using standard tasks. The individual visibility files are 
stacked. Each of these stacks contains typically ~ 300 000 
visibility points, each one averaged over 64 seconds. The 
source subtraction is implemented in the aperture plane us- 
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Figure 2. CLEANed VSA mosaics built up from pointings CrB-A,B,C,D,E,H,I, J,K. In the un-source-subtracted mosaic (left panel) the 
circumferences indicate the primary beams FWHM (2.1°) of the nine pointings, and the crosses the positions of the monitored sources 
with measured fluxes above 50 mjy. The right panel shows the CLEANed source-subtracted mosaic. The synthesised beam FWHM is 
shown in the bottom left corner ((=s 11 arcmin). The noise level is practically uniform across the mosaic at a level m 12 mjy/beam. 
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Figure 4. CLEANed and source-subtracted VSA maps corresponding to pointings CrB-F (left panel) and CrB-G (right panel). The 
nomenclature is the same as explained for Figure 131 



ing AIPS task UVSUB, and then the maps for each individ- 
ual pointing are produced. We used natural weighting, the 
most appropriate in this case since it produces the high- 
est sensitivity and given that the sampling of the aperture 
plane is practically uniform. From the nine central point- 



ings we produced the mosaiced map using standard aips 
tasks. LTESS and STESS were used to make a map in primary 
beam-corrected units, followed by a sensitivity map, which is 
practically uniform around the centre and fails in the outer 
regions. The two mosaics were divided to give the signal- 
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Figure 3. The central region of the CLEANed source-subtracted 
mosaic shown in right panel of Figure[5] We indicate the locations 
of the four negative decrements. We also show the positions of the 
clusters extracted from NED. Five-pointed stars are used for the 
CrB-SC clusters, triangles for the other Abell clusters, squares 
for the Zwicky clusters, and diamonds for the clusters present in 
other catalogues. The centre of the supercluster is indicated with 
a pentagon. Solid contours show positive flux values, while dashed 
contours show negative flux values, where all contours correspond 

to 1.5(7 n . 



to-noise, and the result multiplied by the central value of 
the sensitivity map in order to recover units of mjy/beam. 
In order to avoid loss of signal to noise in the overlapping 
regions, and taking into account that the synthesised beams 
of the nine inner pointings have similar shapes and orienta- 
tions, the CLEAN mechanism (i.e. the deconvolution of the 
synthesised beam) has been applied directly to the mosaic, 
instead of to the individual pointings. To this end, we used 
the synthesised beam corresponding to the central pointing 
CrB-B. We placed CLEAN-boxes around the strongest fea- 
tures, and CLEANed down to a depth of 2<r„. In the case 
of the un-source-subtracted mosaics, CLEAN-boxes were also 
placed around the identified radio sources. The same method 
was applied to the clean map corresponding to pointing 
CrB-F. In the case of the CrB-G map, as the signal-to-noise 
ratio is lower and it is not easy to disentangle the real fea- 
tures from the artefacts, we placed a CLEAN-box around the 
region encompassing the primary beam FWHM, and CLEAN 
was applied down to a depth of 3a„ in this case. 

The CLEANed mosaics, before and after source subtrac- 
tion, are shown in Figure [5] In Figure [3] we present a larger 
scale plot of the regions of interest of the CLEANed and 
source-subtracted mosaic. The maps resulting from the in- 
dividual pointings CrB-F and CrB-G are shown in Figure[I] 
The sensitivity values obtained in each individual pointing 
are shown in Table [5] and the noise level of the overall mo- 
saic is 12 mjy/beam. 



4.1 SZ effect from known clusters in the CrB 
region 

In Table we show the basic data of the eight clusters be- 
longing to the CrB-SC. The X-ray luminosity and electron 
temperature ar e given for the four cluste rs included in the 
BCS catalogue febeling et alJll998t Eoooll : we have also in- 
cluded the more distant clusters A2069 and A2073, for which 
there are also data in BCS. Since we do not have information 
on the parameter (wh ich describes the slope of the den - 
sity profile in a /3-model jCavaliere fc Fusco-Femiand ll976) % ) 
and the c ore radius from X-rays, we have used the relation 
given by iHernandez-Monteaeudo fc Rubino-Martfnl (l200-4h 
in terms of the X-ray luminosity (L x ), 



Art 



-(0.24 ±0.06) 



- (0.1-2.4 keV) 



10 44 h^erg s" 1 



mK 



(2) 



in order to make an estimate of the expected central SZ 
decrement at the VSA frequency 



(AT S z)o = -43.545fj,K 



keV 



.10" 



Mpc 



/[I + (r/r c f 



-3/3/2 



d(r/r c ) 



(3) 



where T e is the electron temperature, n e o the central elec- 
tron density and r c the core radius. 

Making use of the assumption that the inter-cluster 
(IC) medium is a perfect monoatomic gas in thermal equi- 
librium, it follows that the temperature is proportional 
to M a (l + z), with a = |, while the core radius is 
proportional to M~b(1 + z)~^. Hence, we can obtain 



a scaling relation of the type r c 



r c0 M 6 (T e /T e0 )", 
where T e0 = 10 s K, M = 10 15 Afo, and r c0 = 
0.13ft - 1 Mpc iMarkevitch. Forman. Sarazin. fc Vikhlininl 
1998). These have been used to obtain a rough estimate 
of the core radius of each of the clusters for which we have 
X-ray information. To convert to angular sizes, a cosmol- 
ogy Om = 0.25, Q a = 0.75, h = 0.7 3, as derived from 
the last VSA results llRebolo et al 1l2004Tl . was used. Making 
use of these values, and assuming /3 = | , we can simulate 
the response of the VSA to each of these six clusters for 
which we have X-ray information. The brightness temper- 
ature map is converted into flux density, multiplied by the 
VSA primary beam response, and Fourier transformed to 
obtain the simulated aperture-plane response of the VSA to 
the SZ decrement. We must also consider the primary beam 
attenuation for clusters displaced from the pointing position. 
The cluster profiles in the aperture plane for the different 
observing baselines are shown in Figure It is worthwhile 
comparing this plot wi th that presented in the VSA SZ work 
llLancaster et alJl2004h . which shows profiles with a higher 
amplitude, and similar shapes. Since, in that study, over a 
selected sample of rich Abell clusters and with similar noise 
levels, VSA data have shown clear SZ detections in five out 
of the seven observed clusters, we do not expect a high- 
significance detection of SZ effect from the single clusters in 
the CrB-SC region. Moreover, in this case the decrements 
are affected by the primary beam response, except in clus- 
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Figure 5. Predicted SZ flux profiles in the aperture plane for four 
CrB-SC clusters plus other two clusters in the region, as a function 
of the projected baseline length. The thick solid line represents the 
VSA weighting function (arbitrary units), computed from these 
observations. 



ters A2019 and A2124, in which the pointing centres and 
cluster coordinates coincide. 

On the other hand, the simulated visibilities for each 
cluster, computed as explained above, were convolved with 
the synthesised beam of the closest observation in order to 
measure the central flux decrements expected in the maps. 
These values are shown in TableQunder the symbol AS£ rod . 
If we compare the values with the uncertainty in our mea- 
surements, we expect to have SZ decrements from single 
clusters in the region typically with a confidence level be- 
tween 1 and 2a, where a includes now the primordial CMB 
(pcMB ~ 19 mjy/beam), the thermal noise (a n ) and the 
residual sources (a sour ~ 6 mjy/beam) contributions. This 
is exactly what we obtain: comparing the cluster flux values 
from the final maps (see last column in Table , we find 
excellent agreement for those clusters that have a prediction 
for the flux. In all ten clusters analysed, we get negative flux 
values, and in the two most luminous clusters in the CrB- 
SC region, A2061 and A2065, we have 2a detections. These 
two clusters are located close to larger decrements (see Fig- 
ure [3J. Moreover, according to Figure (which includes all 
galaxies lying within ±0.02 of the redshift reported in NED 
for each Abell cluster), A2065 is located in the region with 
the highest galaxy projected density. We can now combine 
all these individual measurements to obtain a statistical de- 
tection — 24.4 ± 9.2 mjy/beam (2.7c). If we also include in 
this weighted mean those values from A2069 and A2073, we 
find -24.8 ±8.1 mjy/beam (3.0a), for all the Abell clusters 
listed in Table [T] 

Finally, another point worth noting is the presence of a 
5a n (signal-to-noise level) decrement (decrement I, see Fig- 
ure^ inside pointing CrB-I, located in a region with a large 
concentration of galaxy clusters, including A2089 and sev- 
eral Zwicky clusters. Their individual SZ effects could have 
an important contribution to its total decrement. 



4.2 The origin of the negative spots 

The most remarkable feature of the mosaic is the existence 
of two prominent negative spots, at signal-to-noise levels 
of 6 and 10, situated inside the primary beam FWHM of 
pointings CrB-B and CrB-H, respectively. Both decrements 
are extended; indeed they cover an area equivalent to ~ 3 
VSA synthesised beams. Coordinates, maximum negative 
flux densities and brightness temperatures of these strong 
decrements are listed in Tabled To convert the flux densities 
into brightness temperatures, we have used the Rayleigh- 
Jeans expression: 



Trj 



A 



-S 



(4) 



2kB^lsb 

where A is the observing frequency, fes is the Boltzmann 
constant, and Q s i, is the solid angle subtended by the syn- 
thesised beam. 

These decrements are located in regions with no known 
clusters of galaxies (see Figure Decrement B is located 
near the centre of the supercluster, and slightly towards 
the north is the more distant Abell c luster A2069, which 
ISmall. Ma. Sargent, fc Hamilton! lll998l) have considered as 
a possible supercluster because of the large number of galax- 
ies it contains. Decrement I and another negative feature 
(decrement J) located towards the east of the mosaic, inside 
pointing CrB-J, have similar flux densities to that of decre- 
ment B, but with a lower significance (signal-to-noise of 5 
in both cases). Note that this latter feature is located in the 
same position as the radio source 1532±2919, so the error 
bar of its measured flux (51 ± 12 mjy, see Table 0, may 
introduce an additional uncertainty. 

We searc hed for possible n egative features in the 
WMAP maps IIBennett et al . 2003) in the region of the most 
intense decrement. Given the angular resolution of the VSA 
extended configuration, we selected the WMAP W-band 
(94 GHz, FWHM=0.21°) map. From this, we can predict 
the observed visibilities as seen by the VSA, enabling us to 
build a "dirty" map. Given the sensitivity of the WMAP 
map (~ 170 /iK per 7 arcmin pixel in this region, the cor- 
responding error would be ~ 108 fiK in a VSA synthesised 
beam of 11 arcmin), we expect to see decrement H only 
marginally at the 2cr n level. In the processed W-band map 
we do indeed see a negative spot in the region of decrement H 
in which the minimum temperature value, ~ — 146± 108 /iK, 
is located at 15 h 21 m 19" ± 28°55'27" (J2000). Note that al- 
though there is an offset of the observed centres of the spots 
between the VSA and WMAP of « 10 arcmin, this is con- 
sistent with the instrumental resolutions. Given the noise 
levels involved, we can not conclude anything about the na- 
ture of the decrement, i.e. we cannot disentangle at la n 
level whether the spectrum between 33 GHz and 94 GHz 
corresponds to primordial CMB (the decrement should have 
temperature —230 /j,K at 94 GHz in this case) or to SZ decre- 
ment (in which case the signal should have a temperature 
-185 fiK at 94 GHz). 

In principle, decrements B and H, given their sizes and 
amplitudes, could be either extraordinarily large primordial 
CMB spots, or SZ signals related to either unknown clusters 
or to diffuse extended warm/hot gas in the supercluster. It 
should be noted that only one third of the galaxies in the 
CrB region are linked to clusters, so SZ contributions could 
be expected in places where there are no catalogued clus- 
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Table 4. Coordinates, flux densities and temperatures of the two largest decrements found in the mosaic. We 
quote both the thermal noise ((J n ) and, within parenthesis, the quadrature of the primordial CMB, thermal noise 
and residual sources + a n + ffour) 2 ) error bars. The last column indicate the percentage over the 5000 

simulations with minima under each decrement flux density (within parenthesis, we also show this value when we 
apply a smoothing to the data; sec Section 4.2.1 for details). 



RA (J2000) 


DEC(J2000) 


AS„ (mjy/beam) 


AT flJ (fiK) 


% below AS„ 


Decrement B 15 25 21.60 
Decrement H 15 22 11.47 


+29 32 40.7 
+28 54 06.2 


-70 + 12 (±24) 
-103 + 10 (±22) 


-157 ±27 (±53) 
-230 ± 23 (±49) 


37.82 (54.16) 
0.33 (0.74) 



ters. As shown in Figures Q and |3 there are few known 
galaxies around decrement H, while decrement B is close 
to a large concentration. This is important since the posi- 
tion of galaxies c ould trace the warm/hot gas distrib ution 
(see, for example. iHernandez-Monteaeudo et al] ]2004l) 'l. In 
the next subsections we explore in detail the three possible 
explanations for the observed decrements. 

4.2.1 CMB anisotropy 

In order to quantify the possible contribution from the pri- 
mordial CMB to these large spots, w e carried out Monte- 
Carlo simulations. Using CMBFAST jSeliak fc Zaldarriaeal 
1996) we generated a CMB power spectrum with a cosmo- 
logical model defined by the following parameters: Qb = 
0.044, Q, M = 0.25, Q.k = 0.75, h = 0.73, r = 0. 14, n s = 0-97 
as de rived from the most recent VSA results <| Rebolo et al] 
l2004l) . plus Tcmb = 2.725 K, = jBennett et al]l2003t) . 
and assuming adiabatic initial conditions. For each decre- 
ment we used this power spectrum to carry out 5000 simu- 
lations of VSA CMB obser vations, follow i ng th e procedure 
explained in Section 3.3 of ISavage et al] ]2004l) . and using 
the aperture plane coverage from each pointings CrB-B and 
CrB-H as templates. Each visibility point contains the CMB 
and thermal noise contributions, plus the confusion level in- 
troduced by the residual radio sources below the subtraction 
threshold of 50 mjy. Thus V = V cm b + V„ + V SOU r, where 

ns 

Vsour(u) = Sie 2lIU ' Xi , (5) 
1 = 1 

where u is the position vector of each visibility point, and Si 
and Xi the flux and the position vector in the map plane of 
the i-th source. The fluxes Si are generated using the source 
count derived from equation yielding a total of ns = 894 
sources distributed between 1 and 50 mjy in a region within 
2° of each pointing. The positions x~l are randomly and uni- 
formly distributed inside this region. 

In the last column of Table 2] we quote the percentage 
of realizations in which the minimum CMB flux value is be- 
low that found in the real map in the two cases in question. 
According to these results the probabilities of decrements B 
and H being caused by primordial CMB (adding the thermal 
noise and the residual sources components) are 37.82% and 
0.33% respectively. Note that here we are considering only 
the intensity: the angular size of the spots is disregarded. In 
order to account for this, we have applied Gaussian smooth- 
ing (with a FWHM approximately equal to the FWHM of 
the spots) to both the simulated and the real maps, and 
in this case the probabilities of the decrements are higher 
(see also Table [IJ. Hence, particularly for decrement H, the 



low probabilities are mainly due to the negative flux density 
values rather than the corresponding angular sizes. 

We calculate the standard deviations of all pixels lo- 
cated within the 2.F-FWHM of the primary beam in the 
5000 realizations, in order to estimate the confusion level 
introduced by the primordial CMB, the thermal noise and 
the residual sources components added in quadrature (<r = 
V a cmb + a n + °iour, where we have a cmb w 19 mjy/beam, 
cr n « 11 mjy/beam and a aour w 6 mjy/beam). From this, 
we find that decrements B and H are deviations at 2.95 and 
4.62o\ These results clearly make a primordial CMB fluc- 
tuation alone an unlikely explanation for decrement H. We 
now focus solely on decrement H. 

Following iRubifio-Martm fc Sunvaevl ]2003t) we per- 
formed a fluctuation analysis. In Figure HJ we present a log- 
arithmic plot of the histograms (i.e. the P(D) functions) of 
pixels values from the 5000 realizations and also using the 
real data from pointing CrB-H. In the case of the realiza- 
tions, the la error bars (containing the primordial CMB, the 
thermal noise and the residual sources) are displayed. The 
plot shows a clear excess in the real data compared with 
the simulations at flux densities below ~ — 30 mjy/beam, 
caused by the presence of decrement H. The excess in the 
positive tail of the curve comes from the two hot spots lo- 
cated towards the north and the east of decrement H in 
the final mosaic (see Fig. |3J. These positive features could 
be either real hot spots in the sky or artefacts due to the 
presence of decrement H itself. (Note also that the response 
of an interferometer has zero mean, so the presence of a 
strong negative feature in the map could enhance the neigh- 
bouring positive spots). It may be that these bright spots 
are enhanced by the presence of radio sources not identified 
in the extrapolation. However, these positive features are 
not present in the map of the region constructed using only 
VSA baselines > 125A, suggesting that these structures are 
extended and not point sources. In order to check if these 
positive spots are due to the sidelobes introduced by decre- 
ment H in the convolution with the synthesised beam, we 
CLEANed the region and find that we obtain lower residu- 
als by placing a CLEAN-box around the negative spot than 
when the CLEAN-box is placed around the positive spots. 
Furthermore, with the box around the negative hole, the in- 
tensity of the positive spot is reduced by ~ 12 %. This is 
clear evidence that these bright spots are enhanced by the 
sidelobes of the synthesised beam. Using the CLEANed map, 
the positive tail of the histogram is strongly reduced, while 
the negative tail is clearly dominant. 

We also have computed the power s pectrum of t he mo - 
saic, following the method described in lScott et al] ]2003j) . 
This is shown in Figure [7] along with that obtained from 
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Figure 7. Power spect rum com puted from the CrB mosaic, in comparison with that obtained from the most recent VSA observations of 
primordial CMB fields iDickinson et al .120041) . The deviation at £ fa 550 is due to the presence of decrements B and H (see Section 4.2.1 
for details). 
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Figure 6. Histograms comparing the distribution of flux den- 
sities inside the 2.1°-FWHM of the primary beam in the 5000 
simulations (thin curve) and in the real data of pointing CrB-H 
(thick curve). We also show the la error bars of the simulations 
(dotted curve). 

the most recent primordia l CMB observations with the 
VSA iDickinson et alJ l2004) for comparison. The compar- 
ison shows a 2.3<r deviation from the pure primordial CMB 
behaviour at I « 550. This is as expected, since this scale 
corresponds to the size of the large decrements. Indeed, we 
find that when we compute the power spectrum using data 
from pointing CrB-H alone, this offset increases to 3. Oct. 
Moreover, as a consistency check, when we remove point- 
ings CrB-B and CrB-H, which contain decrements B and 
H, the resulting power spectrum is then compatible with 
primordial CMB fluctuations on all scales. 

4-2.2 SZ effect from clusters of galaxies 

In order to explore the possible SZ contribut i on, w e fol- 
lowe d the formali s m uti lized by_|Hold«r^^^J (|2000j) (see 
also iKneissl et all <200ll) and iBattve fc Welled J2003T) ') to 

compute the number of clusters in this region potentially 
capable of producing a decrement like H. We estimate the 
mass a cluster would need to have in order to generate an SZ 



effect at least as intense as a given threshold. In this case 
the threshold is the amplitude of decrement H, although, 
to give an idea of the total number of clusters that would 
be detected in the region, we also considered the confu- 
sion level introduced by the quadrature sum of the primor- 
dial CMB, the thermal noise and the residual sources (a w 
23 mjy/beam). We then calculate the number of clusters 
per unit redshift in the Universe with masses above this in a 
solid angle e qual to the entire surveyed region (« 24 deg 2 ), 
by using th e IPress fc Schechterl jl974l) (PS) mass function 
and also the lSheth fc Tor mcn ( 1999J) (ST) mass function, as 
this provides a better fit to the simulations (note however 
that simulati ons also omit much physics). We h ave adopted 
fgas = 0.12 iMohr. Mathiesen. fc Evrardlll999l) for the gas 
mass fraction of the cluster, S c = 1.69 iPe eblcs 1980) for 
the critical overdensity, erg = 0.84 (B ennett et al.ll2003fl for 
the spectrum variance in spheres of radius 8h _1 Mpc, and 
the same values as stated above for the other cosmological 
parameters. 

Figure|5]shows the derived threshold mass and the num- 
ber of clusters per unit redshift versus redshift. By integrat- 
ing this last curve we obtain the total number of expected 
clusters in the surveyed region, which is presented in Tabled 
The number of SZ clusters that should be detected above 
the confusion level introduced by the primordial CMB, the 
thermal noise and the residual sources in the whole surveyed 
region is 3 and 5 respectively for the PS and the ST mass 
functions. They should be located at mean redshifts 0.049 
or 0.091, which is close to the mean redshift of the CrB-SC 
(« 0.07). Note that in our maps we have detected three clus- 
ters in the region above this confusion level (see Table 0. 
On the other hand, the number of clusters that could pro- 
duce a decrement at least as intense as H is respectively only 
0.3 or 0.4 for the PS and the ST prescriptions. Therefore, 
the probability of decrement H of being caused entirely by 
a single cluster is very low. However, if we take into account 
that a decrement such as H due to SZ may be enhanced by 
the primordial CMB, this probability must be higher. Note 
finally that these statistics are valid only for a random patch 
of sky. As we are observing a selected overdense region in 
the direction of a supercluster, the actual probabilities must 
be higher. To obtain a more precise result we may turn to 
N-body simulations, although our analysis is sufficient to get 
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Table 5. Maximum threshold masses, their redshifts, number of detectable clusters in the whole survey, and the most 
likely redshift for clusters producing a SZ effect larger than either the confusion level introduced by the primordial 
CMB, the thermal noise and the residual sources and the minimum flux density of pointing H. We have used both 
the PS and the ST mass functions. 



Threshold 


Maximum 


z where M Hm 


Total number 


z where dN/dz 




M Um (1O 14 M ) 


is maximum 


of clusters 


is maximum 








PS 


ST 


PS ST 


|Min-H|=103 mjy/beam 


27.65 


1.01 


0.3 


0.4 


0.021 0.035 


\j "Ll + + a sour ~ 23 mjy/beam 


11.27 


1.01 


3.2 


5.4 


0.049 0.091 
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Figure 8. Minimum SZ effect detectable mass versus redshift 
(top) and number of expected clusters per redshift unit versus red- 
shift (bottom), using as thresholds either the sum in quadrature 
of the standard deviations of the primordial CMB, the thermal 
noise and the residual sources (thin curve), and the flux density 
of spot H (thick curve). The solid and dashed lines correspond 
respectively to the PS and ST mass functions. 

an order of magnitude estimation, and to show that these 
probabilities are low. 

4-2.3 Diffuse extended warm/hot gas 

As indicated in Section 1, superclusters may be reser- 
voirs of diffuse warm/hot gas, which in principle might 
produce a detectable impr int in the low-energy X- 
ray bands. IZap pacost a et alJ i2002l) . and more recently 
IZappacosta et alTl200^) . have claimed an excess of X-ray 
emission from diffuse structures in a high galactic latitude 
ROSAT field, and the detection of a significant correla- 
tion between the ROSAT-PSPC pointings and the Miinster 



Redshift Survey of galaxies in the Sculptor supercluster re- 
gion. The SZ effect may provide another tool to detect this 
gas, thanks to the long path lengths of the photons across 
these large-scale structures. Using the COBE-DMK data, 
iMolnar fc Birkinshawl il998f) were unable to find evidence 
of an SZ imprint in the region of the Shapley supercluster. 

Since decrement B, according to our studies, can be 
explained through a primordial CMB decrement, we shall 
focus on decrement H, and explore whether it could have 
been built up from a concentration of diffuse warm/hot gas 
(T e < 1 keV). Such a structure producing a detectable SZ 
effect could also leave an imprint on the less energetic X-ray 
bands. We have analysed the ROSAT XRT/PSPC All-Sky 
Survey iSnowden et al.lfl997l) map, corresponding to the R6 
band (0.73-1.56 keV), in order to search for correlated X- 
ray diffuse emission in the region. The R6 map does not 
show an excess of emission in the position of the decrement 
(sec Figure |UJ. The signal in the region of the decrement is 
« 80 — 140 x 10 -6 counts/s/arcmin 2 , which comes mainly 
from the background. In order to perform the ROSAT- VSA 
correlation, we have adopted two different methods. 

Firstly, we have applied a pixel-to-pixel com- 
parison of the ROSAT-K6 map and the real space 
CLEANed mosaic following the metho d descr ibed in 
iHernandez-Monteaeudo fe Rubino-Martml J2004I) . This 
method considers the brightness temperature measured by 
a CMB experiment such as the VSA as the sum of different 
components: a cosmological signal Tcmb, the template we 
want to measure M (in this case, the thermal SZ component 
traced by ROSAT), instrumental noise N, and foreground 
residuals F. The total signal measured at a given position 
on the sky is then modelled as T — Tcmb + otM + N + F, 
where a measures the amplitude of the template induced 
signal. If all the other components have zero mean and 
well known correlations functions, and C denotes the 
correlation matrix of the CMB and noise components, then 
the estimate of a and its statistical error are 



a ~ MC~ 1 M T ' ° a ~ V MC-^MT ■ W 

The derived value of the correlation between ROSAT and 
the VSA mosaic is a = -0.28 ± 0.74 (jlK/X (X denotes 
the units of the R6 map, which are counts s _1 arcmin -2 ), 
which means that there is no significant anticorrelation, as 
would be expected if it were the case that the decrement is 
produced by a diffuse gas distribution. 

We repeated this analysis but in visibility space, by 
predicting the expected visibilities from the ROSAT map 
as seen by the VSA, and performing a similar comparison 
between the observed and predicted values. In this case, we 
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Figure 10. Constraints placed on the overdensity (expressed in terms of baryon density over the mean baryon density) and the electron 
temperature in the centre of the assumed gas concentration (see Section 4.2.3 for details) that generates decrement H. We take into 
account that the X-ray emission must be below that in the ROSAT-R6 map (solid line), and that the depth of this structure should 
be lower than 40 Mpc (dashed line), or lower than 20 Mpc (dotted line). Filled regions indicate the acceptable combinations in the 
parameter space (for 40 Mpc). The factor f represents the fraction of the decrement that is produced by the SZ effect. 
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Figure 9. ROSAT-R6 (0.73-1.56 keV) iSnowden et alJll997ft 

map centred on pointing CrB-H. We indicate the positions of 
all the Abell clusters in the region, emphasizing the CrB-SC clus- 
ters with a star. The units indicated in the upper flux bar are 
10 -3 counts/s/arcmin 2 . The circle indicate the FWHM of the 
primary beam of pointing CrB-H, while the ellipse corresponds 
to the position of decrement H in the VSA mosaic. There is no 
clear excess of emission in this region. 



find a = -0.21 ±0.57 fiK/X for the CrB-H field. The signif- 
icance of this value was derived by computing the dispersion 
of the q value when performing a series of rotations of one 
map with respect to the other prior to the computation of 
the predicted visibilities. 

Assuming the distance of the supercluster (~ 
200ft -1 Mpc, derived from the average redshift of the mem- 
ber clusters (z) = 0.0725), the physical transverse size of 



decrement H is w 2h~ Mpc. At typical WHIM electron tem- 
peratures , i.e. T e = 0.01 keV to T e — 1 keV , a homogeneous 
and spherical gas distribution of this size producing a cen- 
tral SZ decrement of 230 /J.K would have such a high n e that 
it must have collapsed, with subsequent virialisation. This 
should also produce detectable X-ray emission, which is not 
present in the ROSAT-K6 map. In a similar case of a C MB 
decrement without X-ray emission, iKneissl et all lll99Sl) dis- 
cuss various possibilities, but find it difficult to reconcile 
observations with the process of structure formation. A fil- 
ament pointing towards us seems in some ways one of the 
more attractive options. We may therefore consider a larger 
structure with a lower density so that the path length is long 
enough to produce a detectable SZ effect (ATsz oc f n e T e dl) 
without significant X-ray emission (Sx-ray oc J nlT^ 2 dl). 
We can estimate the peak SZ intensity caused by this hypo- 
thetical filamentary structure aligned along the line of sight, 
with a depth L, a central electron density n e , and electron 
temperature T e , at the VSA frequency of 33 GHz: 



AT / = -21.3 \xK 



Mpc 



He 



. 10 _3 cm" 



keV 



(7) 



Here, ATq = 230 /iK is the peak SZ decrement. Given that 
there could also be a contribution from primordial CMB 
anisotropies, we have introduced the factor / which accounts 
for the fraction of the decrement due to the SZ effect. For 
simplicity, we have assumed isothermality and that the elec- 
tron density is homogeneous along the line of sight. We have 
assumed that the depth of this structure must be lower than 
the maximum separation along the line of sight between 
clusters in the core of the CrB-SC, which is 29/i. -1 Mpc 
= 40 Mpc. These assumptions set constraints on the over- 
density (expressed by the baryon density over the universal 
mean baryon density) and T e which are represented in Fig- 
ure 1101 by the dashed (L < 40 Mpc) and the dotted lines 
(L < 20 Mpc). 

We also require the Bremsstrahlung emission to be low 
enough that the signal does not leave a detectable imprint on 
the ROSAT-K6 map. We use equation0to derive a relation 
between L and T e , and integrate along the line of sight. We 
consider a background level of 10 -4 counts/s/arcmin 2 (see 
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Figure [5J as an upper limit for the X-ray emission in the 
region of decrement H. Scaling with the signal of the Coma 
cluster in the same band (63 x 10 -4 counts/s/arcmin 2 ), we 
obtain the constraint represented by the solid line of Fig- 
ure 1101 Filled zones in the parameter space show the ac- 
ceptability regions for structures with lengths of less than 
40 Mpc. These plots show that for / = 1 there is no rea- 
sonable combination of n e and T e to explain the whole 
decrement. Only values of / < 0.7 are able to explain 
the decrement. For instance, a filament with a temperature 
T e « 0.7 keV, a depth L w 40 Mpc, and a baryon over- 
density 5pb/{pb) ~ 850 would produce a X-ray signal of 
fa 0.6 x 10 -4 counts/s/arcmin 2 , which is almost a factor 
2 below the background confusion level of the ROSAT-H6 
map. A structure having these parameters would also pro- 
duce a peak SZ effect of « 115 pK, which is half of the 
amplitude of the decrement (/ = 0.5). Also, a remarkable 
alignment is required to achieve this fraction of brightness 
contrast. 

If instead of the value Spb/{pb) ~ 850, which is com- 
puted for the centre of the filament, we assume an aver- 
age value of Spb/{pb) ~ 450 over the entire structure, and 
a 2 x 2 Mpc 2 square for the transverse shape of the fila- 
ment, the total gas mass enclosed in it would be M gas ~ 
5 x 1O 14 M0. This overdensity appears to be difficult to 
explain in the light of current hydrodynamical galaxy for- 
mation simulations, which predict ovcrdcns ities ~ 10 time s 
lower for these large-scale structures <Cen fc OstrikeJl999D . 
If we assume that the dark matter is well traced by the bary- 
onic matter, this leads to an overdensity of ~ 110 in terms 
of the total matter, in contrast to the typical values for non- 
bound supercluster scales with 5p/p c ~ 5 — 40, for structures 
that have just become virialised with Sp/p c ~ 200, or inside 
the typical Abell clusters with Sp/p c ~ 1000. We note that 
the derived gas mass contained in this hypothetical filament 
is comparable to the total baryonic mass in the known CrB- 
SC clusters, and would represent ~ 10% of the total baryonic 
mass of the supercluster. 



5 CONCLUSIONS 

We have selected the wholly unusual Corona Borealis su- 
percluster region for CMB observations with the VSA ex- 
tended array at 33 GHz. The structures detected in the mo- 
saiced map of the core of the supercluster have significant 
contributions from primordial CMB anisotropics. However, 
these maps show negative flux values at the positions of 
the ten richest clusters in the region, and the two most 
X-ray luminous CrB-SC clusters, A2061 and A2065, each 
produce 2a decrements. When we combine the flux values 
from the positions of the CrB-SC member clusters, we ob- 
tain a statistical SZ detection with a significance of 2.7cr 
(24.4 ± 9.2 mjy/beam). If the clusters A2069 and A2073, 
which have higher redshifts, are added to the sample, we 
obtain 3. Oct (24.8 ± 8.1 mjy/beam). In addition to this, the 
mosaic shows, in a region with a high density of Zwicky 
clusters, a large negative feature, at a signal-to-noise ratio 
of 5, which could be enhanced by the SZ effects generated 
by these single clusters. 

This mosaic shows two strong and extended negative 
features with flux densities —70 ± 12 mjy/beam (—157 ± 



27 pK), and -103 ± 10 mjy/beam (-230 ± 23 pK), 
and positions 15 h 25 m 21.60 s + 29°32'40.7" (J2000) and 
15' l 22 m 11.47 s + 28°54'06.2" (J2000), respectively. These 
decrements are located at positions where there are no 
known clusters. The first one is placed near the reported 
centre of the supercluster. In order to disentangle their ori- 
gins, we have considered the possibility of large primordial 
CMB fluctuations, or SZ signals related to either unknown 
clusters or to diffuse extended warm/hot gas in the inter- 
galactic medium of the supercluster. The primordial con- 
tribution has been explored by performing simulations of 
primordial CMB observations with the VSA. These have 
shown that the size and intensity of decrement B are consis- 
tent with primordial CMB fluctuations, whereas in the case 
of decrement H the probability is only 0.33%. Decrement H 
is marginally detected in the W-band map of the WMAP's 
first-year dataset, but the current sensitivity does not allow 
us to disentangle the spectral behaviour of this decrement. 

In order to investigate the second possibility we have 
predicted the number of clusters capable of producing a 
decrement such as H. To this end we used the Press-Schecter 
(PS) and also the improved Sheth-Tormen (ST) prescrip- 
tions to describe the density of collapsed objects in the Uni- 
verse. This study showed that the number of random clusters 
in the VSA observed fields massive enough to produce this 
feature is 0.3 and 0.4 respectively for the PS and the ST 
formalisms. Note that this probability has been computed 
assuming that the whole decrement is due to the SZ effect. If 
we take into account that there could be a primordial CMB 
contribution, this probability becomes higher. 

The ROSAT-K6 data show no evidence of X-ray emis- 
sion in the region of decrement H. This fact sets constraints 
on the electron density and temperature that could have a 
hypothetical intercluster warm/hot gas distribution capable 
of producing a decrement as deep as H. If it were in a small 
volume, then the electron density would be high enough to 
produce detectable X-ray emission in the typical ranges of 
WHIM temperatures (0.01 — 1 keV). Hence, we hypothe- 
sise that the decrement could be caused by a large filament 
aligned in the direction of our line of sight. We have ex- 
plored the possibility that only some fraction of the spot is 
due to an SZ effect, the primordial CMB being responsible 
for the rest. In this case, a filament with a length below 
the size of the supercluster, i.e. < 40 Mpc, a temperature of 
0.5 — 0.8 keV, and a baryon density in its centre between 500 
and 900 times the mean baryon density in the local Universe, 
produces a SZ effect close to one half the central flux den- 
sity of the decrement, and its X-ray emission would be low 
and obscured by the background in the ROSAT-K6 map. It 
would contain a gas mass of ~ 5 x 1O 14 M0, which is com- 
parable to the total baryonic mass contained in the CrB-SC 
member clusters. If the value 3 x 10 16 Mq is assumed for the 
total mass of the supercluster, this filament would hold the 
~ 10% of the total expected baryonic mass of the superclus- 
ter. However, the required overdensity of such structure is 
~ 10 times higher than the predictions from N-body galaxy 
formation simulations for these large-scale filaments. But we 
stress that: i) we are observing at a very atypical region, and 
ii) N-body simulations need more physics. 

In summary, we are confident that our measurements 
do show excess decrement, and to explain decrement H we 
require a combination of primordial CMB fluctuations with 
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either an SZ effect from an unknown cluster or from a large- 
scale filamentary structure, which would hold a significant 
fraction of the total baryonic mass of the supercluster. It 
is worthwhile to carry out multi-frequency observations of 
decrement H in order to disentangle these possibilities. 
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